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ABSTRACT 

We have performed extensive two-dimensional magnetohydrodynamic simulations to study the am- 
plification of magnetic fields when a supernova blast wave propagates into a turbulent interstellar 
plasma. The blast wave is driven by injecting high pressure in the simulation domain. The in- 
terstellar magnetic field can be amplified by two different processes, occurring in different regions. 
One is facilitated by the fiuid vorticity generated by the "rippled" shock front interacting with the 
background turbul ence. The resulting turbul ent fiow keeps amplifying the magnetic field, consistent 



with earlier work (Giacalone & Jokipii 2007). The other process is facilitated by the growth of the 
Ray leigh- Taylor instability at the contact discontinuity between the eject a and the shocked medium. 
This can efficiently amplify the magnetic field and tends to produce the highest magnetic field. We 
investigate the dependence of the amplification on numerical parameters such as grid-cell size and 
on various physical parameters. We show the magnetic field has a characteristic radial profile that 
the downstream magnetic field gets progressively stronger away from the shock. This is because the 
downstream magnetic field needs a finite time to reach the efficient amplification, and will get further 
amplified in the Ray leigh- Taylor region. In our simulation we do not observe a systematic strong mag- 
netic field within a small distance to the shock. This indicates that if the magnetic-field amplification 
in supernova remnants indeed occurs near the shock front, other processes such as three-dimensional 
instabilities, plasma kinetics and/or cosmic ray effect may need to be considered to explain the strong 
magnetic field in supernova remnants. 

Subject headings: shock waves - magnetic field - turbulence 



1. INTRODUCTION 

Powerful shocks associated with supernova remnants 
(hereafter SNRs) sweeping through the interstellar 
medium (ISM) are remarkable high-energy phenomena 
in astrophysics. It is widely believed that these high- 
Mach number shocks are the sources of galactic cosmic 
rays with energies up to at least 10^^ eV. SNRs are also 
sources of strong radio and/or X-ray emissions. In these 
high-energy processes, the magnetic field is of great im- 
portance. Moreover, it provides information on energetic 
charged particles, which are presumably accelerated by 
the supernova shocks. 

The ISM is known to be turbulent. Measurements 
of the ISM radio-wave scintillation have established the 
existence of large-scale density turbulence which has a 
Kolmogorov-like power spectrum spanning more than ten 
decades of s patial scale with an outer scale of several 
parsecs (e.g., Lee fc Jokip ii"1976^ 'Armstrong et al. "19811 
' Rickett||l99Ql lArmstrong'et al . 1995; Minter fc Spangler] 
1996D. This has been called "the big power law in the 



sky." (see, Spangler 2007) The galactic magnetic field 
is observed to be a lew micro-Gauss and has uniform 
and fiuctuating c ompone nts that are rough ly in equipar- 
tition (e.g.. Beck et al .||7996; Mint er fc Spangler] |199"6 ; 
|Han et al .l2Q04[ ). The turbulent magnetic field can mter- 
act with tne shock waves, distorting their surfaces lead- 
ing to shock ripples ( Neugebauer fc Giacalone|2005 [) and 
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the enhanced downstream magnetic fiuctuations (Lu et 
al. 2009 |. It is also important for efficient particle accel^ 
eration ( Giacalone""2005j Jokipii fc Giacalone 2007 Guo 
fc Giacalone 2010). Turb ulence in the upstrea m medium 
has also been considered ([Balsara et a l.|[2 00l| to explain 



the irregular and patchy emission niorpho logy observed 
in SNRs (e.g., [Anderson fc Rudnick||1996D . 

Recently, it has been interred Irom observations that 
the magnetic field in young SNRs is strongly enhanced 
to a magnitude much greater than the compression given 
by the shock jump condition. For example, by assuming 
that the so- called X-ray "thin rims" seen in several young 
SNRs (e.g., Bamba et al.|200 5^ are caused by shock accel- 
erated electrons rapidly losing energy in strong magnetic 
field through synchrotron radiation, the associated mag- 
netic field may be more than 50 — 100/iG in o rder to 
explain the thickness of the "thin r i ms'' (Berezhko et al. 



2003 



i] lyolk et al.|[2QQ5l lBallet|[2QQ6| [PiZizot et al.| |2006[ ) 

BNK shocks with higher shock speeds propagatim 



I^br SNK shocks with higher shock speeds propagating 
in more inhomogenous media such as Cas A and Tycho, 
the downstream magnetic fields are inferred to be more 
enhanced. "Thin rim s" are also seen in radio emissions 
( Reynoso et al.||1997 ), which cannot be explained by the 
electrons losing energy in strong magnetic field. This 
indicates some other mechanism, e.g., d ecay of magneti c 
fiuctuations may need to be considered (Pohl et al.|2005|). 
Further downstream of the shock, the magnetic field may 
possibl y be even higher than the region right behind the 
shock ( |Vink fc Laming||2003 ). The morphology of X-ray 
emission in SNRs shows filamentary structure and rapid 
time variation, which indicates th at the mag netic field 
could be as high as ImG ( |Uchiyama et al. 2007), over two 
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orders of magnitude higher than the background fluid. It 
should be pointed out this rapid variation of synchrotron 
emission could be well reproduced in the case of strong 
magnetic fluctuations ( |Bykov et al.||2QQ8 ). Since the ef- 
fect of turbulence is not fully understood at this point, 
the actually magnetic-field amplification factors in young 
SNRs remain uncertain. 



Beh fc Lucek| ( |2QQ1 ) and |Bell| ( |2QQ4[ ) proposed that 
cosmic rays acceleratea by the S^R forward shock waves 
provide a current that leads to an instability that can 
amplify the magnetic field close to the shock front. Nu- 
merical simulations show the evidence of this instability, 
however, it is found to easily saturate and the a mplifica- 



tion f actor may b e limited (e.^., [ Ri quelme & Spitkovsky 
' Jo kipiij (|2QQ7) proposed 



20091). Recently, Giacalone 

an alternative mechanism, in which thelnteraction be- 
tween the warped shock front and large-scale density fluc- 
tuations produces fluid vorticity downstream of strong 
shocks. That fluid vorticity can stretch, distort and 
amplify the magnetic field. The magnetic-field ampli- 
fication in this mechanism relies on the dynamics of 
a magnetized fluid rather than the cosmic- ray kinetic 
physics. It is interesting to note, however, I Balsara "et 
[aL (2001| performed three-dimensional MHD blast wave 
simulation with moderate resolution which does not show 
magnetic-field enhancement larger than 50/iG, whereas 
two-dimensional Cartesian geometry simulations with 
high numerical reso lution give strong amplification (Gi- 



[acalone fc Jokipii|2007| jlnoue e t al. 2009) with maximum 
values larger than 100/iG. This discrepancy warrants fur- 
ther investigation. 

In addition, a critical constraint from several observed 
"thin rims" is that the required magnetic-field amplifi- 
cation should occur within a narrow distan ce < 0.01 pc 
of the su pernova shocks (for a review, see iReynolds et 



al.||2011 ). This is supported by some coincidences of X- 
ray "thin rims" and shock locat ions inferred from Ha 
observations ( Winkler et aH2003 ) and radio polarimetry 
( Gotthelf et al.||2001| )l 'I'his places an important con- 
straint on various held amplification models, although it 
should be noted that there is no one-to-one correspon- 
dence between X-ray "thin rims" and inferred shock lo- 
cations. To our knowledge, this constraint has not been 
taken into account previously in comparing models and 
simulations with observations. 

In this study, we perform a series of two-dimensional 
ideal MHD simulations with high spatial resolution to 
study strong supernova blast shock waves propagating 
into the ISM containing pre-specified large-scale density 
and magnetic fluctuations. We investigate how the am- 
plification depends on a variety of parameters, including 
the explosion energy, the level of background turbulence, 
and the numerical resolution, etc. The paper is orga- 
nized as follows. In Section 2, we describe our numerical 
model and simulation setup. We present the simulation 
results in Section 3, and in Section 4 we summarize and 
discuss our results. 



2. BASIC CONSIDERATION AND NUMERICAL 
MODEL 

We have performed a series of two-dimensional ideal 
MHD simulations with high-order numerical schemes. 



The ideal MHD equation can be expressed as: 

S,p + V-(pu) = (1) 

dtipu) + V • {pun + (P + ^)I - ^) = (2) 

S,(^) + V-((^+p+|^)u-^^^) = (3) 

S^B- V-(uB-Bu) = (4) 

where p is the plasma density, u is the fluid velocity, B 
is the magnetic field, p is the gas pressure, and E is the 
total energy density, which is defined as: 
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We use a central finite- volume scheme on overlapping 
grid cells ( Liu et al.,^200 8) to solve the ideal MHD equa- 
tions. In particular, a high-order divergence-free re- 
construction for the magnetic fiel d that uses the face- 
centered values has been employed ( Li [2008 , 2010). The 
magnetic field is advanced with a high-order constrained 
transport (CT) scheme to preserve the divergence- free 
condition to machine round-off error. The overlapping 
cells are natural to be used to calculate the electric-field 
flux without a spatial- averaging procedure, and hence we 
can achieve the higher-order accuracy (3rd-order or 4th- 
order). The central schemes do not need time-consuming 
characteristic decompositions, and are easy to code and 
be combined with un-split discretization of the source 
and parabolic terms. The overlapping cell representation 
of the solutions is also used to develop more compact re- 
construction and less dissipative schemes. For solutions 
that contain discontinuities, e.g., shocks and/or contact 
discontinuities, we apply a non-oscillatory hierarchical 
reconstruction (HR) to remove the spurious oscillations 
and achieve high resolution near the disco ntinuities. The 
HR limitine; we have used ( |Li |2008| 



2010) requires infor- 



mation only from the nearest neighbor cells and it does 
not require characteristic decomposition. The numerical 
dissipation introduced by the HR limiting is enough to 
damp out all the artificial oscillations near the shocks and 
no extra artificial viscosity is needed. To further improve 
the computational efficiency and reduce the numerical 
dissipation for the smooth flow, we develop a shock de- 
tector to flag the cells near the shock and perform HR 
only to those cells. The d etails of the wh ole algorithm 



have been documented in (Li 20081 2010). Our method 



has been verifled to achieve the expected order of accu- 
racy and have very low numerical dissipation. The high- 
order, low-dissipation, and divergence-free properties of 
this method make it an ideal tool for MHD turbulence 
simulations. 

We model the simulation in a two-dimensional Carte- 
sian coordinate (x, y) with uniform grids. The size of the 
simulation domain corresponds to x Ly = 30pc x 30pc 
in all the simulations. A supernova blast wave is driven 
by the initial injection of thermal pressure and mass in 
a small circular region at the center of the simulation 
box. The physical parameters used in the simulations 
are listed in Table 1. To calculate the injection energy, 
we must set a volume for the injection region. We as- 
sume a cylinder with radius of r = OApc and length in z 
direction is Lz = 0.8pc, which gives a volume of 0.402pc^. 
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Initially, the cylinder region has a high pressure pin^ and 
a total mass about 3 solar masses. For most cases con- 
sidered here, we assume that the injected internal energy 
in the central region is 1.5 x 10^^ erg. A grid number of 
X fly = 4000 X 4000 is typically used for the simula- 
tion domain. This initial s etup for blast waves i s similar 
to the simulation made by Balsara et al. (2001), except 
that we use two-dimensional simulations with higher res- 
olutions to study the magnetic-field evolution in young 
SNRs. The density and magnetic field in the background 
ISM consist of an average component and a turbulent 
component. We take the interstellar background plasma 
mean number density to be no = 1 cm~^. The average 
magnetic field Bq is along the y direction, whose magni- 
tude is given in Table 1 for different cases. In our model, 
we assume a constant initial background ISM pressure po. 
Its corresponding temperature Tq for each case is listed 
in Table 1. For the fiuctuating components, we assume 
that both magnetic field and density fiuctuations have a 
two-dimensional Kolmogorov-like power spectrum: 



P (X 



1 + {kLc)^/^ 



(6) 
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This is consistent with observations of interstellar tur- 
bulence ( [Arm strong et al. 1995 Chepurnov & Lazarian 
I. The turbulence is generated by summing a l arge 



number of discrete wave modes with random phases ( Gi- 
acalone fc Jokipi'i||1999 ). In all the simulation cases the 
coherence length of the background turbulence is chosen 
as Lc = 3 pc. The random component of magnetic field 
is given by 



SB{x,y) = ^ ^y CB27rknAknPB{kn){s'i^ OnX - COS 6ny) 
n=l 

X exp(z COS OnknX + i sin OnknH + ^0ni7) 

where PB{kn) is the power for wave mode with 
wavenumber /c^, random propagation angle —1 < 
cosOn < 1 and random phase < 0n < 27r. is a 
constant used to normalize the wave amplitude. In this 
study, the total fiuctuating magnetic field is taken to be 
SB'^ = Bq in ah cases. 

The density fi uctuations satisfy a log-normal probabil- 
ity d istribution (Burlaga & Lazarus 2000 Giacalone & 
Jokiph,2007J : 



n{x,y) = noexp(/o + (5/) 



(8) 



in which /o is a constant used to determine the mean 
density and 5f is given by a similar expression as the 
turbulent part of the magnetic field 

^/(^, y) = En^i y^Cf2nknAknPf{kn) 

X exp(i cos OnknX + i sin OnknV + ^^n) (9) 

In Runs 1 - 4 we consider the effect of the amplitude 
of upstream density fiuctuations. In Runs 5 and 6 we 
examine the effect of different value of ^o- In Runs 7 
and 8 we examine the effect of different temperatures of 
the ISM. We have also simulated two different explosion 
energies in Run 9 and Run 10. In Runs 11 and 12, we 
examine the effect of the simulation grid-cell size. Runs 



1 and 2 are the representative runs with and without 
pre-specified background turbulence. 

3. SIMULATION RESULTS 

We now describe the simulation results on the interac- 
tion between the supernova blast wave and the turbulent 
ISM, along with the magnetic-field evolution downstream 
of the supernova shock. We will first present the results 
from our representative Run 1. Then we will discuss the 
spatial distribution of magnetic-field amplification, fol- 
lowed by the effects of numerical resolution and other 
effects such as the background turbulence and explosion 
energy. 

3.1. Two Regions of Magnetic- field Amplification 

Figure [l] shows the time evolution of (a) the total ki- 
netic energy summed over the whole simulation domain, 
(b) the average radius of the blast wave, and (c) the av- 
erage speed of the shock for Run 1. The average radius 
is estimated using the area of the high pressure region 
for each snapshot of pressure in the simulation, and the 
average speed of the shock is calculated by taking run- 
ning differences between the averaged radii. In the be- 
ginning of the simulation, the region with high density 
and high pressure expands and drives a shock propagat- 
ing into the turbulent medium. The kinetic energy in- 
creases sharply during the expansion and reaches about 
1.15 X 10^^ erg, i.e., ~ 77% of the injected explosion en- 
ergy is converted to kinetic energy at its peak. After that 
the swept-up ISM slows down the eject a so the kinetic 
energy decreases slowly. In all the cases we have sim- 
ulated, the total energy is conserved within a degree of 
10~^ during the simulation time. The shock speed de- 
creases from about 8700 km/s to about 3000 km/s at the 
end of the simulation. The corresponding Alfven Mach 
number changes from ~ 850 to ~ 290. The radius of the 
remnant roughly follows r oc t^'^, which is consistent with 
the self-similar solution fC hevalie'r||1982 ). The size, age 
and shock speed are roughly consistent with the obser- 
vations of young SNRs, though there is a wide spread in 
these quantities observationally. After a few times of the 
ejected mass being swept by the supernova blast wave, 
the shock speed is expected to slow down and settle into 
an "Sedov" phase in which r (x.t^-^ typically after a few 
thousand years. Since in three-dimension the interstel- 
lar mass swept by supernova shock is proportional to 
rather than in our two-dimensional simulation, the 
shock speed in later times (after several thousand years) 
is likely overestimated. In this work we mainly focus on 
the evolution of young SNRs. 

Figure |2] (a) and (h) display the total magnetic-field 
energy ana the maximum magnetic-field strength in the 
simulation box for Run 1. It is shown that the magnetic- 
field energy increases to a level much larger than the 
initial magnetic-field energy. At the end of the simula- 
tion, the magnetic-field energy reaches to ~ 0.47 x 10^^ 
erg. These results are qualitatively consistent with pre- 
vious studies using the planar shock waves ( Giacalo ne fc| 
Jokipii|2007]|Inoue et al.|2009| . The maximum magnetic 
held can rapidly increase to 250/iG within 500 years and 
finally goes up to about 550/iG in 3000 years, much larger 
than what is expected from the jump condition. But the 
detailed analysis shows that the spatial location of the 
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maximum magnetic field is typically not near the imme- 
diate downstream region of the shock. We will discuss 
this result in more detail later. 

Figure [3 shows the snapshots of (a) the velocity magni- 
tude, (b) the magnitude of magnetic field, (c) density and 
(d) temperature at t = 1600 years for Run 1. It is shown 
that the velocity field of the blast wave is highly irregular. 
The shock surface is rippled as regio ns with different den- 



sities pass through the shock front (Giacalone & Jokipii 
2007[) . The flow at the rippled shock transition produces 



strong transverse and rotational flow downstream of the 
shock wave. It can be seen from Figure [s] (b) that the 
magnetic fleld downstream of the blast wave is strongly 
amplifled. We flnd the ampliflcation is c losely related 
to the down stream vorticity production (^Giacalone &' 
Jokipii|2007|). This flow patten stretches and distorts the 
field lines of force, which leads to a small-scale dynamo 
process. 

In addition, different from the shock amplification, we 
find that the magnetic field in the interface region be- 
tween the eject a and the shocked medium is also strongly 
enhanced by the Rayleigh- Ta ylor instability (RTI) a t 
the contact discontinuity (e.g., I J un fc Norman| |1996a|bl). 
It appears that the magnetic field can be ennanced to 
~ 300/iG in the region where RTI is important. Figure 

t(c) shows "fingers" of enhanced density resulting from 
e RTI in the interface region, which is quite different 
from the mechanism discussed above, which occurs just 
behind the shock. In fact, the magnetic-field amplifica- 
tion by the RTI process contributes almost equally to 
the total enhanced magnetic energy and the maximum 
field strength as shown in Figure |2] The nonlinear de- 
velopment of RTI stretches the magnetic field and causes 
strong amplification. This process has been studie d ex- 
tensively (e.g., see the earlier num erical studies by Jun et| 
al.|1995t|Jun fc Normanjl996a|b[ ). The RTI amplification 
process can be studied m more detail when a more real- 
istic initial eject a profile are taken into account. Though 
beyond the scope of this paper, it will be worthwhile to 
see if there are possible observational signatures in the 
RTI-amplification region. 

3.2. Spatial Dependence of the Magnetic-field 
Amplification 

As discussed in the Introduction, some observed X- 
ray "thin rims" and their coincidence with the inferred 
shock locations have suggested that magnetic field is am- 
plified at or within a short distance to the supernova 
remnant shock front. It is thus imperative to investi- 
gate whether the MHD simulations can reproduce this 
feature, although to directly compare to X-ray observa- 
tion would also require the simultaneous computation of 
cosmic-ray variation, which is not included in this study. 

In Figure [4] we present the results at t = 600 years 
for Run 12, which has the highest available resolution. 
It shows the probability distribution functions (PDFs) 
of magnetic field, taken within a thin region of 0.15 pc, 
0.3 pc and 0.45 pc behind the shock front, respectively. 
In addition, the PDF of the downstream magnetic field 
within a distance of 0.3 pc of the shock front for Run 2 
(without background turbulence) is also plotted for com- 
parison. The PDFs presented here and in other figures of 
this paper are normalized by the respective volume from 
which the distribution is taken. It appears that the mag- 



nitude of magnetic field increases with the distance from 
the shock front. Specifically, for the rim within 0.15 pc 
downstream, the maximum magnetic field reaches only 
70/iG. Also the region with magnetic field higher than 
30/iG only occupies about 0.8% of the rim. 

To explore further the amplification of the magnetic 
field in both the shock downstream region and the 
Ray leigh- Taylor region, we plot the average magnitude 
of magnetic field as a function of radial distance from 
the shock front at different times from 200 years to 2000 
years for Run 12, which is shown in Figure |5[ Together 
with the spatial distribution of magnetic field at these 
times (not shown here) that is similar to Figure |3j we 
can separate approximately the two different magnetic- 
field amplification regions. In Figure [5j the right side 
of the dashed lines is typically dominated by the shock 
amplification and the left is by the Ray leigh- Taylor in- 
stability. It is shown that Ray leigh- Taylor amplification 
produces much larger magnetic field than those from the 
shock amplification. Detailed analysis shows that, lo- 
cally, the amplified magnetic field in the Rayleigh- Taylor 
region can easily reach several hundred micro-gauss. 

It can be seen from the downstream magnetic-field evo- 
lution that the amplification of magnetic field has not 
reached a saturation. This is probably not surprising, 
given the fast transit time of the SNR shock and the rel- 
atively young age of the SNRs. The gradual increase of 
the magnetic-field amplitude further downstream from 
the shock is consistent with the picture that turbulence 
has longer time to amplify the fields. Before the magnetic 
field reaches Rayleigh- Taylor region, its maximum value 
is only about 200 — 300 micro-Gauss, much less than the 
amplified magnetic field in Rayleigh- Taylor regions. 

3.3. Effects of Numerical Resolution 

Since the field amplification at the shock front is closely 
related to the vorticity generation, the numerical resolu- 
tion is expected to play an important role. Numerical 
modeling of thisprocess must resolve the vortical mo- 
tions. In Figure |6] (a) we plot the normalized probability 
distribution function (PDF) of the magnitude of mag- 
netic field downstream within a distance of 0.3 pc of the 
shock front for three different grid resolutions. Run 1 and 
2 have 4000 x 4000, Run 11 and 12 have 2000 x 2000, and 
8000 X 8000, respectively. Figure |6](b) shows the evolu- 
tion of the total magnetic-field energy for these cases. 
It can be seen that higher resolution runs give higher 
total magnetic-field energy and, even for the highest res- 
olution of Run 12, the total amount of magnetic energy 
in the simulation domain has not converged. More de- 
tailed analysis shows that the kinetic energy and mag- 
netic energy have not reached equal partition in small 
scales. Note that the total magnetic energy at the end 
of the simulation (~ 3 x 10^ yrs) is still a small fraction 
(< 10~^) of the injected explosion energy or the available 
kinetic energy in the simulation. 

In Figure [6] (a), by plotting the PDFs within a small 
downstream region close to the shock front, we can get a 
more clear view of the effects of the numerical resolution 
in the shock amplification process near the shock front. 
The green curve from Run 2 represents the shock ampli- 
fication without the background ISM density turbulence, 
whereas Run 1, 11, and 12 represent the shock amplifica- 
tion being significantly enhanced when the background 
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ISM turbulence is present. The purple curve represents 
the initial, un-shocked background ISM magnetic-field 
distribution in the same region as shown in the cases 
that include ISM density turbulence. The higher resolu- 
tion Run 12 obviously produces a greater volume frac- 
tion with higher magnetic field and a larger maximum 
magnetic field in the downstream region of the shock, al- 
though one might argue that the difference among these 
runs is not very large, at least for the region within 0.3 
pc of the shock at this particular time. 

In Figure [g] (b), however, the difference among differ- 
ent resolutions seems to be much bigger. Such differ- 
ences are mostly caused by the region where RTI ampli- 
fies the magnetic field (see Fig. [s] (b) for field distribu- 
tion), because the total magnetic energy is a summation 
of both the shock-amplified region and the R TI refflon. 
This is consistent with some previous studies ( Jun et al. 
|T995) where numerical resolutions are shown to be im- 
portant as well. Such requirements for high resolution 
might also explain the previous fin ding that mametic- 
field amplification is not as strong (Balsara et al. 2001) 
when it is difficult to employ very high numerical resolu- 
tion in three-dimensional MHD simulations. We have 
also checked the long ter m evolution (up to 3 x 10^ 
years) of SNRs similar to Balsara et al. (2001) using 
two-dimensional high-resolution simulation and we con- 
sistently find strong magnetic field larger than 100/iG in 
the downstream region. 

Figure [6] demonstrates that both the shock and RTI 
amplification of the magnetic field depends on the nu- 
merical resolution (perhaps more strongly for the RTI re- 
gion). Figure [t] further supports this conclusion where we 
plot the vorticity distribution of the shocked flows for the 
low resolution (Run 11) and the high resolution (Run 12), 
at a time that is the same as in Figure |3] It is not surpris- 
ing to see that more small scale structures are developed 
with much larger vorticity magnitudes in Run 12. Note 
that the two magnetic-field amplification regions can be 
approximately separated spatially (at least at this time). 
The small scale structures in vorticity (signifying turbu- 
lence) are developed in the immediate downstream of the 
forward shock region whereas both relatively large (i.e., 
the density "fingers") and small scale vorticity features 
are produced in the RTI region. 

3.4. Other Effects 

Figure [S] shows the effects of different background tur- 
bulence amplitudes. In Figure [sl (a), the PDFs of mag- 
netic field within 0.3 pc behind the shock front for Runs 
1 - 4 at t = 600 years are shown, respectively. Figure 
M (b) represents the evolution of magnetic-field energy 
lor these cases. It is seen that the larger amplitude den- 
sity fiuctuation tends to lead to stronger magnetic-field 
amplification. For Run 2, it can be seen that the effect 
of RTI strongly enhance the magnetic field. The total 
magnetic-field energy can even exceed the case of Run 4 
{^Srfi /tIq = 0.3). Detailed analysis shows the growth of 
RTI is somewhat suppressed by turbulence which makes 
the magnetic-field energy in Run 4 less than Run 2 at 
late times. 

We also examine the effect of the magnitude of initial 
interstellar magnetic field in Run 5 and Run 6. Compar- 
ing these two cases with Run 1, we find the shock ampli- 



fied magnetic field is nearly proportional to the magni- 
tude of initial magnetic field. This is consistent with the 
fact that the amplified fields have not reached saturation. 
The effect of the temperature of background medium is 
also examined in Run 7 and Run 8. We find that the 
different temperatures do not yield any strong difference 
in the downstream magnetic-field evolution. 

Figure [9] is similar to Figure |8l but for different initial 
explosion energies, £^ = 1.5 x 1(F^ erg, 4.5 x 10^^ erg, and 
4.5 X 10^^ erg, respectively. The time frames are chosen 
so the shock radii are roughly the same, the magnetic- 
field distribution before amplification is therefore roughly 
the same. We can see the magnetic-field amplification is 
stronger for the case of higher explosion energy. This is 
because the central region drives a stronger shock which 
generates stronger vorticity in the downstream region. 
The magnetic energy evolution in these three cases fol- 
lows the same trend. 

4. DISCUSSION AND CONCLUSION 

The inferred strong magnetic field in young SNRs is 
a significant result and can be important in the high- 
energy processes including particle acceleration and ther- 
mal/nonthermal emissions. The origin of this process, 
however, is still under debate. In this work we study 
the interaction between a supernova blast wave with a 
turbulent upstream medium which contains density and 
magnetic-field fiuct nations. The vorticity produced at 
the rippled shock front can stretch and distort the mag- 
netic field lines, and this leads to a strong magnetic-field 
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amplification do wnstream ( Giacalone fc Jokipii[2007 
oue et al. 2009). Using two-dimensional MHD siniula^ 



tions of a blast wave, we confirm the key features of this 
process. Based on our simulations, we conclude that the 
increase of magnetic field is dependent on shock speed 
and background density turbulence amplitudes. Further- 
more, the numerical resolution used in the simulations 



can play an imp ortant role as well. Previous work (Bal- 
sara et al. 2001) using three-dimensional MHD sirnula- 
tion with moderate resolution shows no magnetic-field 
amplification beyond 50/iG. Here we show the magnetic 
evolution downstream is sensitive to the resolutions used 
in the simulation. For high resolutions, the simulations 
allow rapid growth at small scales, this leads to efficient 
field amplification. 

Furthermore, we find that there are two different pro- 
cesses and spatial regions where magnetic field is ampli- 
fied. One is associated with the shock amplification im- 
mediately downstream and the other is associated with 
the RTI at the interface between the eject a and the 
shocked medium. 

However, in our simulations, we did not observe a sys- 
tematic strong magnetic field within a thin region im- 
mediate downstream of the supernova shock. For ex- 
ample, using the results of the highest resolution case, 
within 0.15 pc downstream of supernova shock, we ob- 
serve only about 0.8% region which has magnetic field 
larger than 30/iG. This lack of strong magnetic field can 
be understood as the downstream dynamo process re- 
quires an efficient stretching to produce strong magnetic 
field. The time scale for the growth of magnetic field 
depends on the eddy turnover time. Only after a certain 
time can the field get sufficient amplification. If the thin 
rims (0.01 — 0.1 pc) observed in young SNRs are indeed 



6 



caused by the electrons losing energy in strong fields (~ 
several hundred /iG), some other processes such as three- 
dimensional instabilities, plasma kinetics or the effect of 
cosmic rays might be needed to explain the magnetic- 
field amplification in young SNRs. 

We note that in observation there is no one-to-one 
correspondence between X-ray "thin rims" and inferred 
shock locations. In fact many "thin rims" and filamen- 
tary structures seen in X-ray observation can hardly be 
related to shock front due to observation limitation. Fur- 
ther understanding about the relationship between the 
small scale X-ray structure and shock locations is needed 
to further constraint and distinguish the different mech- 
anisms for amplification of magnetic field. 

We also note that the two-dimensional simulation in 
this study could be significantly different from three- 
dimensional simulation. It is known that in two- 
dimensional simulation the dynamics of MHD flow is 
very different from that for three-dimensional simula- 
tion. For example the inverse cascade of enstrophy can 
causes strong intermittency due to two-dimensional ef- 



fect (Biskamp 2003). There are several recent studies 
show that magn etic field amplification behind high-mach 
number shock fln oue et al .''2011a* b|) and in Rayleigh- 
Taylor region (Stone ^ Gardiner 20^ can still operate 
in three-dimensions, which confirm tne results found in 
two-dimensional simulations. Also, in three-dimensional 
simulation the MHD flow could develop other types of 
instabilities with larger growth rates. Further three- 
dimensional MHD simulation with high resolution will 
be useful in confirming the conclusions of this paper. 

ACKNOWLEDGEMENT 

This work was supported by the LDRD and IGPP 
programs at LANL and by DOE office of science via 
CMSO. Computations were performed using the insti- 
tutional computing resources at LANL. JG, JRJ, and 
FG also acknowledge partial support from NASA grant 
NNX10AF24G. JRJ is also supported in part by NASA 
grant NNX11AB45G. We acknowledge useful discussions 
with Dr. Hao Xu and Dr. Federico Fraschetti. 



REFERENCES 



Anderson, M. C, &; Rudnick, L. 1996, ApJ, 456, 234 
Armstrong, J. W., Cordes, J. M., & Rickett, B. J. 1981, Nature, 
291, 561 

Armstrong, J. W., Rickett, B. J., & Spangler, S. R. 1995, ApJ, 
443, 209 

Ballet, J. 2006, Advances in Space Research, 37, 1902 

Balsara, D., Benjamin, R. A., &; Cox, D. P. 2001, ApJ, 563, 800 

Bamba, A., Yamazaki, R., Yoshida, T., Terasawa, T., & Koyama, 

K. 2005, ApJ, 621, 793 
Beck, R., Brandenburg, A., Moss, D., Shukurov, A., &; Sokoloff, 

D. 1996, ARA&A, 34, 155 
Bell, A. R., &; Lucek, S. G. 2001, MNRAS, 321, 433 
Bell, A. R. 2004, MNRAS, 353, 550 

Berezhko, E. C, Ksenofontov, L. T., Volk, H. J. 2003, A&A, 412, 
Lll 

Biskamp, D. 2003, Magnetohydrodynamic Turbulence, Cambridge 
University Press 

Brouillette, M. 2002, Annual Review of Fluid Mechanics, 34, 445 
Burlaga, L. F., & Lazarus, A. J. 2000, J. Geophys. Res., 105, 2357 
Bykov, A. M., Uvarov, Y. A., &; Ellison, D. C. 2008, ApJ, 689, 
L133 

Chepurnov, A., & Lazarian, A. 2010, ApJ, 710, 853 
Chevalier, R. A. 1982, ApJ, 258, 790 
Giacalone, J., & Jokipii, J. R. 1999, ApJ, 520, 204 
Giacalone, J. 2005, ApJ, 624, 765 
Giacalone, J., &; Jokipii, J. R. 2007, ApJ, 663, L41 
Giacalone, J., & Neugebauer, M. 2008, ApJ, 673, 629 
Gotthelf, E. v., Koralesky, B., Rudnick, L., et al. 2001, ApJ, 552, 
L39 

Guo, F., & Giacalone, J. 2010, ApJ, 715, 406 

Han, J. L., Ferriere, K., & Manchester, R. N. 2004, ApJ, 610, 820 
Inoue, T., Yamazaki, R., & Inutsuka, S.-i. 2009, ApJ, 695, 825 
Inoue, T., Asano, K., &; loka, K. 2011, ApJ, 734, 77 



Ino ue, T., Y amazaki, R., Inutsuka, S.-i., &; Fukui, Y. 2011, 

arXiv:1106.3080 
Jokipii, J. R., & Giacalone, J. 2007, ApJ, 660, 336 
Jun, B.-L, Norman, M. L., & Stone, J. M. 1995, ApJ, 453, 332 
Jun, B.-L, & Norman, M. L. 1996, ApJ, 465, 800 
Jun, B.-L, & Norman, M. L. 1996, ApJ, 472, 245 
Lee, L. C, & Jokipii, J. R. 1976, ApJ, 206, 735 
Li, S., 2008, J. Comput. Phys, 227, 7368-7393 
Li, S., 2010, J. Comput. Phys., 229, 7893-7910 
Liu Y., Shu C.-W., Tadmor E., and Zhang M., 2007, Commun. 

Comput. Phys., 2, 933-963 
Lu, Q., Hu, Q., & Zank, G. P. 2009, ApJ, 706, 687 
Minter, A. H., & Spangler, S. R. 1996, ApJ, 458, 194 
Neugebauer, M., &; Giacalone, J. 2005, Journal of Geophysical 

Research (Space Physics), 110, 12106 
Parizot, E., Marcowith, A., Ballet, J., &; Gallant, Y. A. 2006, 

A&A, 453, 387 
Pohl, M., Yan, H., & Lazarian, A. 2005, ApJ, 626, LlOl 
Reynolds, S. P., Gaensler, B. M., & Bocchino, F. 2011, 

Space Sci. Rev., 46 

Reynoso, E. M., Moffett, D. A., Goss, W. M., Dubner, G. M., 
Dickel, J. R., Reynolds, S. P., &; Giacani, E. B. 1997, ApJ, 491, 
816 

Rickett, B. J. 1990, ARA&A, 28, 561 

Riquelme, M. A., &; Spitkovsky, A. 2009, ApJ, 694, 626 

Stone, J. M., &; Gardiner, T. 2007, ApJ, 671, 1726 

Spangler, S. 2007, Bulletin of the American Astronomical Society, 

38, 162. Presentation online: 

phobos.physics.uiowa.edu/~srs/AAS_07.ppt 
Uchiyama, Y., Aharonian, F. A., Tanaka, T., Takahashi, T., &; 

Maeda, Y. 2007, Nature, 449, 576 
Vink, J., & Laming, J. M. 2003, ApJ, 584, 758 
Volk, H. J., Berezhko, E. G., & Ksenofontov, L. T. 2005, A&A, 

433, 229 

Winkler, P. F., Gupta, G., &; Long, K. S. 2003, ApJ, 585, 324 



Shock Turbulence Interaction 



7 



Run 


n 


V ' u 


^0 






Grids 


Ba7;e(rim) 


^macc (rim) 


\ 


1.0 


0.45 


3 


1.0 X 10^ 


1.5 


4qqq2 


8.7 


48.0 


2 


1 n 


n n 
u.u 


Q 

o 


l.U A lU 


1 ^ 
1 .o 


zLnnn2 
^uuu 


o.u 


23 2 


Q 

o 


1 n 


n 7 


Q 
O 


1.0 X 10^ 








ux .o 


4 


1 n 


u.ou 


Q 
O 


1.0 X 10^ 






o.o 




5 


1.0 


0.45 


]^ 


1.0 X 10^ 


1.5 


4000^ 


2.9 


16.1 


6 


1.0 


0.45 


9 


1.0 X 10^ 


1.5 


4000^ 


25.9 


143.4 


7 


1.0 


0.45 


3 


1.0 X 10^ 


1.5 


4000^ 


8.7 


48.1 


8 


1.0 


0.45 


3 


1.0 X 10^ 


1.5 


4000^ 


8.7 


48.2 


9 


1.0 


0.45 


3 


1.0 X 10^ 


4.5 


4000^ 


9.0 


63.2 


10 


1.0 


0.45 


3 


1.0 X 10^ 


0.45 


4000^ 


8.3 


37.7 


11 


1.0 


0.45 


3 


1.0 X 10^ 


1.5 


20002 


8.3 


35.6 


12 


1.0 


0.45 


3 


1.0 X 10"^ 


1.5 


8000^ 


9.1 


72.8 



TABLE 1 

The parameters used in the simulations. Number density n is in cm~^, magnetic field B is in /iG, temperature Tq is in K, and 

EXPLOSION ENERGY E IS IN ERGS. THE LAST TWO COLUMNS LIST THE AVERAGE AND MAXIMUM MAGNETIC-FIELD STRENGTH IN UNITES OF jlG 
WITHIN A THIN REGION (0.15 PC) DOWNSTREAM OF SUPERNOVA SHOCK AT 600 YEARS. IN ALL THE CASES, THE FLUCTUATING COMPONENT 

OF THE MAGNETIC FIELD IS TAKEN TO BE 53"^ = 
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Fig. 1. — The time evolution of (a) the total kinetic energy summed over the simulation domain, (b) the average radius, and (c) the 
average speed of the shock for Run 1. The initial injected thermal energy is 1.5 X 10^-*^ ergs. 
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Fig. 2. — The time evolution of (a) the total magnetic-field energy and (b) the maximum magnetic-field strength in the simulation 
domain -Bmax, for Run 1. 
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Fig. 3. — Spatial distribution of (a) the velocity magnitude, (b) the magnetic-field strength, (c) density, and (d) temperature for Run 1 
at t = 1600 years. The background ISM turbulence causes the shock front to become rippled. In addition, the interface between the ejecta 
and the shocked ISM experiences the Ray leigh- Taylor-like instabilities. 
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Fig. 4. — The probability distribution function of magnitude of magnetic field downstream with a distance within 0.15 pc (black solid 
line), 0.3 pc (blue dashed line) and 0.45 pc (red dot dashed line) behind the shock front for Run 12 at t = 600 years. The probability 
distribution functions for magnetic field with a distance within 0.3 pc behind the shock front in Run 2 at t = 600 years and the initial 
background magnetic field are also plotted for comparison. 
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Fig. 5. — The average magnetic-field magnitude as a function of radial distance to the shock front from Run 12 at different time from 200 
years to 2000 years. The dashed lines roughly separate magnetic-field regions dominated by the shock amplification and by Ray leigh- Taylor 
convective fiow. 




Fig. 6. — (a) The probability distribution function of magnitude of magnetic field downstream of a distance within 0.3 pc to the shock 
front at t = 600 years for runs with different resolution. The PDF of the background magnetic field in the same region taken from Run 1 
at t = is also plotted, (b) The comparison of total magnetic-field energy evolution for the same runs. 
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Fig. 7. — The color-coded images of vorticity (top) and the magnetic-field amplitude (bottom) for Run 11 (low resolution 2000 x 2000, 
left panels) and Run 12 (high resolution 8000 X 8000, right panels) at t = 1600 years. 
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Fig. 8. — Comparison for different background turbulence amplitudes, (a) The probability distribution function of magnitude of magnetic 
field downstream within a distance of 0.3 pc behind the shock front at t = 600 years, (b) The comparison of the total magnetic-field energy 
evolution for these runs. 
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Fig. 9. — Comparison for different initial explosion energies, (a) The probability distribution function of magnitude of magnetic field 
downstream with a distance within 0.3 pc behind the shock front when the shock radii are roughly the same for these cases, (b) Comparison 
of the total magnetic-field energy for these runs. 



